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Abstract 

We use a covariant supermultiplet theory to determine the primary coupling 
constant associated with several types of two-body meson decay. Despite the 
diverse range of decays considered the primary coupling constant is surprisingly 
uniform. We envisage the extension of the techniques to heavy quark cases, 
including as preliminary examples the calculation of the D* + and D*° total 
decay widths with results 57.7 ± 1.5 KeV and 42.5 ± 2.6 KeV respectively, as 
well as some predictions about D* and B* radiative decays. 
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1 Introduction 



In these heady days of heavy quark effective theories (HQET) we feel it is appropri- 
ate to review a supermultiplet scheme developed in the mid-1960's which has great 
similarities to the HQET in respect of the heavy quark and the accompanying 'brown 
muck'. In doing so we hope to test the covariant supermultiplet theory for light as 
well as heavy degrees of freedom, assessing the extent of symmetry breaking and how 
it manifests itself. We may then apply our techniques to heavy quark examples in 
future work with considerable confidence. 

To determine the matrix elements of currents between hadrons requires knowledge 
of the hadronic wave function in terms of the quark and gluon constituents. A full 
relativistic treatment of such constituents is impossible because of the infinite degrees 
of freedom associated with the quarks and gluons. However, the success of the non- 
relativistic quark model spawned several workers |1|, [| |3|] to construct relativistic 
spinor fields describing pointlike mesons incorporating the correct spin, parity, flavour 
and colour degrees of freedom. Despite later evidence that the mesons were not 
pointlike objects, these group theoretical approaches can be shown to be equivalent 
to the non-relativistic weak binding limit [ffl. 

The basis of our work is a relativistic meson supermultiplet field ||. The wave- 
function describing the meson is dynamically equivalent to a system of two quarks, 
both of which are on-shell and moving at the same velocity. This differs little from the 
heavy quark picture which assumes that a meson with a quark much heavier than its 
light antiparticle partner will have the heavy component almost on-shell and moving 
at the same velocity as the meson because the light 'brown muck' must move with 
the same velocity. This assumption, along with the hope that weak interactions will 
not affect the motion of the heavy quark at small recoil, led to the now popular heavy 
quark symmetry and decoupling (see J/J |8| for a review). 

The reasons for using the supermultiplet scheme are several. Firstly, the scheme 
automatically incorporates the Zweig rules and duality diagrams, so one can easily 
determine the Zweig allowed strong decays. At the same time it incorporates isospin 
and field mixing factors so that one can readily normalise the coupling constant of 
various decays; indeed this makes supermultiplet theory very predictive as one need 
only know a single coupling constant to predict widths of many seemingly unre- 
lated processes. Secondly, radiative decay modes may be examined by combining the 
supermultiplet scheme with the vector meson dominance model. This permits the 
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theory to make some quite accurate predictions about photon mediated decays such 
as uj — > 7r + 7r~. Thirdly, excited mesonic states may be constructed in terms of the 
supermultiplet field |§ so we are able to further broaden its applications. Finally, the 



scheme is easily extended to include the c |T(J and b quark flavoured mesons so that 
we may venture into the heavy quark arena with little modification. 

The recent ACCMOR |JTTJ] and CLEO Collaborations [12[] have renewed interest 
in the D* decays due to two major findings. The D* + total width was measured 
with an upper bound of 131 KeV (signficantly lower than the 1992 upper bound 
of 1.1 MeV), and the D* + — > D + ~f branching fraction appears significantly smaller 
than earlier measurements. Both these findings are consistent with constituent quark 
model |13| and HQET predictions |14]]. We provide similar calculations within the 
supermultiplet framework and reproduce these findings with considerably far ease. 



2 Supermultiplet Field 

Derivation of the wave function |l| proceeds by applying a relativistic boost to the 
rest frame spinor <f) b a {p), odd under parity. This leads to a relativistic spinor <f)^(p), 
satisfying Bargmann-Wigner equations, namely 

4> B A (P) = €M = (*+ m)[ 7 ^ 6 a - 750 5 a02m (1) 

where a, b are flavour indices, a, (3 are spin indices and 5 corresponds to the pseu- 
doscalar nonet and <fi^ to the vector nonet (with p^<p^ = 0). One can show |4| that 
such a relativistic spinor is equivalent to describing the meson as a quark-antiquark 
pair both of which are moving at the same velocity as the meson and are therefore 
both on-shell. 

We use the simplest effective interaction Lagrangian 

C mt = G $f (pO [$g(p2) $£( P3 ) + (P2 <- Ps)) (2) 

as proposed by [|T| to describe the three point coupling between the mesons involved in 
two-body decays. G is a normalization factor and we note that it has the dimensions 



of mass. Such an interaction Lagrangian corresponds to a duality diagram |fLq , |16 
as shown in Figure 1 with three mesons meeting at a vertex (p^ incoming). Here 
flavour labels have been included to show how the flavour is automatically conserved, 
flavour being carried by the line. Contravariant spinor indices correspond to odd 
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parity as they represent the antiquark, while covariant indices have even parity since 
they represent the quark flavour. 

Upon substitution of the supermultiplet field ([I]) into the interaction Lagrangian 
(0), and using the Dirac trace algebra along with the momentum conditions 

Pi + P2 + Ps = 0, pf = m- , = 4>{pi) 

we reduce the interaction Lagrangian to 

= g Y „ fa - PsT <^(1)[0 5 (2),0 5 (3)]> (3) 
+ 9 VVP ^ KX Pi,P2X < M (1)R(2),0 5 (3)} > (4) 
+ 9vw[(P2 - PzYg v °m x + (p 3 - pi) v g a »m 2 + ( Pl - p 2 ) a g^m 3 
+2(p 2 - PzY{pz - Pifipx - p 2 y/ {mi +m 2 + m 3 )] 
<^(l)[^(2),0 ff (3)]> (5) 

where <> stands for a trace over the internal symmetry indices, corresponding to 
a joining of quark lines in a duality diagram. For instance, such a trace for flavour 
indices would expand as 

0; a (l)(0 5fc c (2)0 5c a (3) - 5b c (3)«2)) 

for the vector-pseudoscalar-pseudoscalar (VPP) vertex. The interaction Lagrangian 
contains three distinct coupling constants as at this stage we are alert to the pos- 
siblility that symmetry breaking may affect each piece differently, that is we have 
introduced three separate constants depending on the type of decay. However, full 
supermultiplet symmetry would mean the constants are related in the following way: 

2g VPP = m x g VVP = m x g vvv (6) 

One should also note that the three pseudoscalar meson vertex is not present in the 
interaction Lagrangian, as expected by parity conservation. (This follows automati- 
cally in the supermultiplet scheme because the trace over three 5 fields disappears 
from Equation [2]. ) 

For the moment we only consider strong meson decay so that tree-level calculations 
in perturbation theory given by (§, ^, |^) suffice. Also we apply the general formula 
for a two-body decay, 





where Si is the spin of the parent meson and \(m 2 ,m 2,1713) = m\ + m\ + m\ — 
2m\m\ — 2m\m\ — 2m 2 m 2 . 

Upon substitution of the interaction Lagrangian in the form of Equations (HUH) 
into the decay rate formula (0), we derive the following widths for the various decays: 

IWp = \^\mlmlml)g 2 vpp /A^m\ (8) 
Yy^yp = X 3/2 (mlmlml)g 2 VVp /967rml (9) 
Tv_ V y = \ 3/2 (ml,ml,ml)gl vv y(m u m 2 ,m 3 )/l92nmlm 2 2 ml 1 (10) 

where 

y(m 1 ,m 2 ,m 3 ) = 9 Y (m + mj) 2 (m^ - mj) A - 

\l<i<j<3 

3 _ \ 3 

98 Y ™\- 16 Yl (mi + rrij) 3 +142j]m l 2 . 

i=l l<«<i<3 / i=l 

We now have an adequate formalism for describing various strong interaction de- 
cays amongst ground state mesons. To extend the applications of the supermultiplet 
theory we invoke the ideas of the vector meson dominance model in order to account 
for various electromagnetic interactions of our mesons. To do so we make the usual 
assumption that the coupling between a vector meson flavour singlet and a photon is 
of the form 

g vl {k 2 = 0) = em 2 v /g vpp , (11) 

as shown in Figure 2. When extrapolating away from k 2 = we expect the coupling 
to decrease and as such have denoted the coupling by em v /g' vpp to allow for such 
change. That is, we anticipate g' pp {k 2 ) will vary with k 2 due to intermediate virtual 
particle contributions and its value at k = equals g VVP in Equation ||. 

The vector meson dominance model, used in conjunction with our decay rate 
formulae (PfljlCiD give the following rates for the various processes: 

IW = (m 2 v -Am 2 y/ 2 (l-m 2 /m 2 v ^ 2 (e 2 /g' vpp ) 2 /127r (12) 

= (m 2 v -m 2 P ) 3 (eg vvp /g' vpp ) 2 /96nm v (13) 

IW 7 = (m 2 P -m 2 v ) 3 (eg vvp /g' vpp ) 2 /32irm P (14) 

T P ^ 77 = m P (e 2 g vvp /g' g' ) 2 /64n (15) 



thereby greatly extending the original scope of the supermultiplet scheme. In going 
from our purely strong interaction decay rates to the radiative ones, we have used the 
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gauge invariance of our interaction Lagrangian and simply substituted a mass of zero 
for those vectors connecting with the photon. However, the three vector interaction 
(HD is only gauge invariant for the case 772.2 = so strictly we should only apply 
it to radiative examples for which the virtual vector meson satisfies this condition. 
Unfortunately, since the photon only couples to flavour singlet states the condition 
m 2 = m 3 also implies the daughter vector mesons are identical. Due to the F-type 
coupling between daughter states in the interaction Lagrangian fl5|) such decay widths 
will automatically go to zero. It is for this reason we have not included a^-> V'y 
term above, despite experimental evidence for such (eg. T^^/Y^^w < 2%; although 
our zero width prediction does not conflict with this). We now go one to apply the 
formalism to the ground state mesons. 



3 Supermultiplet Method 



In the standard way we take the pseudoscalar nonet as: 



/ 7T° I Tig I T?Q 

y/2 ^ s/S ^ V3 

K- 



7T 



'V2^ V6^ V3 



K° 



K + 
K° 

V6^ V3 / 



(16) 



where 



7/8 = cos Op 7/ — sin 9p rj 
7/o = sin Op rj + cos Op rj 1 



(17) 
(18) 



as defined in |L7j and Op is the pseudoscalar mixing angle. The vector nonet is 
similarly given by 

p+ K* + \ 



V 



where 



P 

K*~ 



u 8 



'y/2 + VE + V3 



A) 



K 



*o 



K 

V6 ^ V3 I 



(19) 



cos Oy <p — sin Oy u 

Sin Oy (f) + COS Oy Ul 



(20) 
(21) 



and Oy is the vector mixing angle as determined using the Gell-Mann-Okubo (GMO) 
mass relation. We have chosen to go into some detail as our formalism is not the 
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common one adopted by recent publications [IE, [T!|. Conversely, 

= cos 9y uj% + sin 9y ujq (22) 

and the uj field as 

uj = — sin 6 V uj 8 + cos 9 V uj , (23) 



where ujg and ujq masses are determined by the GMO relation ||17|| . The vector mixing 
angle is obtained from 

2(( m| - m 2 )(m 2 - m l)) 1 / 2 
2ml — mi — m 2 



tan2fl y = i 2 8A 2 8 24 



where ?>m\ = 4m 2 K * — m 2 p . 

From the vector nonet ( |i~9|) the octet and singlet fields are expressed in terms of 
the supermulitplet vector as 

^^8 = 0j + <f>* - 20/ 3 , V3 uj = M i + </>* + 0/ 3 . 

Substituting these into Equation ^2| yields 

a/6 = (cos9 v + V^smOv)^^ + <j>£) + (-2 cos fly + v^sinfly)^. 

In the case of "ideal mixing" = M | so that 

cos 6y + a/2 sin 9y = or tanfly = — 1/V2 



leaving us two options for 9y] either —ir/2 < 9y < or 7r/2 < 9y < 71. The first case 
implies cos fly = J2/3, sin fly = so that = —0^1 while the second gives 

the desired result of = |. Thus a suitable solution to Equation Ellis in the range 



>3 

7r/2 < 9 V < n. More generally, the solution to (|24]) is 

. /2((m| - TO 2)( m 2 2))1/2X 

2Q V = tan" 1 U _ 8A „ 8 + nvr 

V 2m s - m i - m i J 

where n is any integer. 

The above arguments for the determination of the vector mixing angle can be 
applied to the pseudoscalar nonet with the substitutions uj$ — > 17$, uj — > 770, — > 
77, uj — > t]'. Using the condition tt/2 < 9y < tt and a similarly derived expression for 
the pseudoscalar angle, — n/2 < 9 P < tt/2, we obtain the equally likely results 

9y = 129.4°, 140.6° 
9 P = -10.5°, 10.5°. 
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The mixing angles we have obtained may seem accurate, but the GMO relation is 
extremely sensitive to an extra small SU(3) symmetry breaking associated with the 
27 representation; a small 27 addition can produce a major modification of the angle. 

With the correct structure now in place, it is a relatively simple process to test 
the supermultiplet theory. We wish to calculate the standard coupling constants 
g VPP and g VVP , examine how similar they are for each process and finally compare 
the supermultiplet prediction (§) of the relation between them. In practice we take 
the decay width and particle masses as input [pl| and determine the coupling constant 
associated with the decay via (||,[D and (p^-[15|). The simplicity of the supermultiplet 
method is that isospin and mixing factors are automatically accounted for. One 
simply chooses an appropriate decay, determines the flavour indices a, b and c using 



matrices (p!6|JT9| ) and then use these in the correct part of the interaction Lagrangian 
([|. [|, or |5|) to determine the normalization factors which arise. For example, in the 
decay p + — > n + it , a = 1, b = 2, c = 1, 2 and upon substitution of the fields into 
Equation [3] one finds gp+Tr+^o = y/2 g VPP so that the coupling constant we determine 
for this decay should be divided by the factor y/2 to obtain the standard coupling 
constant g VVP - This procedure is repeated for all appropriate physical decays. Mixing 

to replace the ideal fields 



is easily accommodated by using the relations ( IT| , |lB , |2T) , 
by the real mesons in the interaction Lagrangian. 

In radiative decays of the type V — > 1 1 we allow for the coupling of the photon to 
the quark. Using the following electromagnetic charge projectors 

/ 2/3 \ 

Ql= 





-1/3 







-1/3 J 



(25) 



we may likewise extract the relevant standard coupling. 

Radiative modes such as V — > P7 require some delicacy in normalising the cou- 
pling constant. To elicit a clear understanding of the method we include an example 
of the procedure for the decay p° — ► 777. Firstly, we recognise p° is a combination 
of </>„i(l) and <#„o(l), so we require terms in the interaction Lagrangian (01) with 



1, b 



1 and a 



a = 1, = 1 ana a = 2, b = 2. For each of these cases we determine the third flavour 
index c such that 4> v (2) is a flavour singlet that couple to a photon (thus c = 1,2). 
Substituting these values into formula f| we pick out the uncharged parts; 



C mt oc g VVP [ 0^(1) {6,i(2),05i(3)} + 0J(1) {0^(2),0 52 2 (3)} 
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„ [P°(1)/Q/P Q (2) g; 8 (2)wq 8 (3) , ^q(3)n 

9vvp I V2 ^' V2 + ^ n/6 ~T V3 } 

-p°(l) / o / -p (2) , m,(2)w»»(3) , 20(3)^1 

*g vvp p°(l)p°(2) [ % (3) + v / 2r/ (3)] 
;cos^p + v^sin^p) p°(l) p°(2) 77(3), 



where in particular we have used relations ( |T7i[l8l) to arrive at the final result. The 
form shows that the coupling between two p° mesons and a pseudoscalar 77 is related 
to the standard VVP coupling by 



9p°p°v = V 1 ( cos Qp + V2 sin dp) g vvp (26) 

The virtual vector meson is immediately identifiable as p°(2), and we must necessarily 
allow for the coupling between this and the photon. Since p° = (uu + dd) / y/2 then 
9p°-y — <?v7/v2 which in turn implies 5^0 pp = V^g' VPP from (P"T|). Subsequently, the 
coupling between a p°, r\ and photon is related to our standard couplings by 

9 P °m = e 9 P a P r,l g'pQpp 

e J| (cos 6 P + v^sintfp)^^ 



V%g' 

v 3 VPP 

4=f (cosfl P + v^sin^p) e ' J] 



V3 g' VPP 

In other decays, it is possible that the radiative mode may proceed via more than 
one virtual vector meson. The above method is still used to determine each virtual 
vector meson contribution and the appropriate linear combination is taken. 



4 Results 

The results are presented in Tables 1 and 2. For clarity these tables include the SU(3) 
factors which we have used to normalize the coupling constant. 

Table 1 summarises the results of our investigation into the coupling between 
a vector meson and two pseudoscalar mesons. The first half of Table 1 displays 
purely strong interaction decays, while the second lists the coupling constant g' 
obtained from vector meson dominance extrapolation. Two important features are 
apparent. Firstly, we have found that the coupling is far more regular than previously 
believed by those persons who deprecate light quark symmetry. Secondly, the form 
of the symmetry breaking is now very clear. As the mass of the parent vector meson 
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increases (as we go down each half of the table), so does the coupling, apparently 

1 /2 

following the simple rule g VPP ~ 0.154 m,y (for my in units MeV). Similarly, the 
mass-shell constants g' VPP follow such a relation, except the constant of proportionality 
is approximately (0.136 ± 0.003)MeV _1//2 by a weighted mean method (and an error 
scale factor of 4; following the Particle Data Group's handling of errors). This result 
complies with the known scaling law behaviour for fy as my — ► oo ||. Thus if the 
matrix elements of quark currents between a given vector meson and the vacuum 
state is defined by 

< | <? 2 7/ t <?i I V >= f v m v e^ 

then it is well known |2(J that fy oc| \ /my 2 as my — > oo. Translating to our 
terminology fy = emyj g' p , we verify this prediction and importantly we find the 
result is also supported in the light meson sector. Admittedly, the ui — *> e + e~ has a 
very high g' , but since the width of uj —>■ fi + fi~ is only known to an upper bound 
(providing a lower bound estimate of g' vpp ) the anomaly remains unsubstantiated. 

Table 2 predominantly lists the results from studying radiative decays to obtain 
estimates of g VVP using vector meson dominance. The first entry in the table is for the 
decay <p — > pn and leads to a direct determination g VVP , not via a radiative transition. 
In fact, we use it to test the supermultiplet prediction 2g vpp = mig VVP , the results 
of which are shown in Figure 3. This figure shows the sum vri\g VVP — 2g VPP plotted 
against vector mixing angle 9y and clearly demonstrates the supermultilpet rule is 
satisfied at 9y « 140.3°, very close to the accepted value 9y = 140.6°, and it is for 
this reason we have used this value in all our calculations. 

In the case of radiative decays, where we know the coupling is related to the 

1 /2 

ratio g VVP / g' vpp , we have used the relationship g' p ~ 0.136 m y , which is well 
supported by the data in Table 1. Importantly, this relation applies to the virtual 
vector meson so that for decays mediated via the ideal field ug we have to use its 
mass of approximately 931 MeV. The data shows that once again the coupling is 
quite regular, but now the symmetry breaking appears to obey a power law relation 
g VVP oc m{ n where 1/2 < n < 3/2. 

5 Predictions 

With a clearer understanding of the effects of symmetry breaking on the coupling 
constant, we may now confidently determine the decay rates for non-Zweig allowed 
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decays. In particular we study the decays u — > it + it~ and (p — > 7r + 7r~, both of which 
are mediated by a virtual photon coupling between the parent vector meson and a 
p meson (electromagnetic mixing). Thus in the case uj — > it + it~ we have the overall 
coupling of 

gun* = e 2 m 2 p g pn7T /g' Wpp {ml - m 2 p )g' ppp 
and for <fi — > 7r + 7r~ we have 

94m* = e 2 m 2 p g pnw /g' (f>PP (ml - m 2 p )g' ppp . 



Using 



g pn7T = (0.1537 ±0.002)mJ /2 

g' ppp = V2 (0.136 ±0.003)my 2 

9'u, PP = (0.136 ±0.003)my 2 / sin 9 V 

9<t> PP = ^ (0-136 ±0.003)mi /2 / cos 9 V 



we predict 

T^ 7T+7T - = (1.66 ± 0.16) x 10~ 2 MeV 
r^ +7r - = (5.88 ± 0.55) x 10" 4 MeV 

which compare favourably with the presently accepted values 

T^ 7T+7T - = (1.86 ± 0.25) x 10~ 2 MeV 
T^ n+n - = (3.5 ± 2.8) x 10" 4 MeV. 

The symmetry breaking effects we have observed also lead to a measurable conse- 
quence in the radiative decays of heavy mesons. We begin by re-examining the decays 
A'* ± — > K ±f y and K*° — > K ^. Experimentally, the K* branching fraction is 

r K*°^K°j/^ K*+^K+j = 2.31 i 0.29 

and allowing for phase space factors this translates into a coupling constant ratio of 

\gK*°K°-y/9K*+K+y\= 1.514 ± 0.095, 

and as such is far from the exact SU(3) ratio of 2. Under the supermultiplet scheme, 
one can show the decays proceed via two intermediate vector mesons, p° and cj 8 . 
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Following the procedure we described for determining the normalisation factors, one 
finds 

(g K *+pOK+ . 9K*+Lu a K+\ / „x 

9k*+k+^ = e — + — (27) 

V 9 p opp 9u> s pp J 

( 9K*°pOKQ gK*+u 8 KQ \ 

9K*°K°-y = e — ; 1 - l . (28) 

V 9 p opp 9u, s pp J 

If one assumes q' is constant then 



"l/yg | -1/V 6 N 
g VPl /3 



/ ±/ V z -IV" / ,/ \ 

fe*+^+ 7 = e — =- + ^ }{9vvp/9 vpp ) 



, -1/V2 -i/VE \ , 
g K *o KOl = e \ — — + I {g VVP /g v 




and we arrive at the exact SU(3) prediction. If instead we use a symmetry breaking 
5^ pp we must substitute g' p op P = y/^Cm 1 ! 2 and g' w&PP = \J~QCm}J 2 in Equations 
(g7D and (H). Thus 

9K-°K0y = _ m pQ 1/2 + m ^ 1/2 / 3 = _j g7 
#if*+K+ 7 m; 1/2 - m~ 8 1/2 /3 

and notice the result is independent of C, the constant of proportionality between 

1 /2 

g' vpp and . Although not matching the experimental result, it is an improvement 
on exact SU(3). Actually, the most satisfactory explanation of the symmetry breaking 
mechanism comes from |^TJ. They attribute the deviation from exact SU(3) to the 
constituent mass difference between the strange and non-strange quarks in the loop of 
a quark triangle diagram. As K* —>■ Kj excite both strange and non-strange quarks, 
such a difference must be accounted for. With these corrections, the experimental 
ratio is found to match theoretical estimates very well. We intend to apply the method 
to heavier meson cases in future work |22| . 



Let us continue to use the "mass variation principle" of g'„ in the heavy meson 
sector. Upon application to the D* and B* mesons we obtain 

Qrr," 1 / 2 -4- m-V2 _L_ Am" 1 / 2 

iS^L = 3m V'% 2 J/ % ~-«> (29) 
g D *+D+~f -3m p0 ' +m a ; s / +4mj^ 



g B -°B° 7 = -3m; 1/2 + m-JJ 2 + 2m- 1 / 2 
gB*+B+ 1 3m~ 1/2 + m^ 8 1/2 + 2rriv 1/2 



-0.34 
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which are significantly different from the exact SU(5) predictions of 



gD*°D°j/gD*+D+-y = 4 

and as such require better experimental data to test the results. 

In addition to these relative decay rate predictions, the supermultiplet scheme 
can be easily adapted to decay width calculations. Scattered amongst Tables 1 and 
2 are various constants determined by extending the supermultiplets to include the 
charm and bottom quark mesons. In particular, using the upper bound of 131 KeV 
for the D* + decay width |11| we have found g VPP < 10. Conversely, we can use our 



knowledge of the effects of symmetry breaking to predict the VPP coupling constant 
for D*+. We find 

g VPP (D* + ) « (0.1537 ± 0.002)(2010) 1/2 = 6.89 ± 0.09 (30) 



surprisingly similar to a heavy quark prediction of 7 ± 1 by ||23|| . We can use the 



g VPP value for D* + to calculate the total decay width of the D* + — > PP channels. 
Using the supermultiplet method we can predict all the possible decays of D* + into 
two pseudoscalars; however, phase space restricts the processes to D* + — ► D tt + and 
D* + — ► D + 7T° so that the width must be 

T D *+^ PP = r D *+^ D + 7r o + T D *+_ iD o 7r + 

= g 2 VPP [^ /2 (rn 2 Dt+ ,m 2 DO ,ml + ) + A^m 2 ^, m^ + , m 2 o)]/487rm^ + 
= 57.7 ±1.5 KeV. 

We compare this with the radiative width D* + — > D +/ y in the following branching 
fraction: 



D* 



>D+7 (g VPl /g VPP ) 2 (m 2 Dt+ - m 2 D+ ) 3 /96nm 3 Dt , 



T D ,+^p P [A 3 / 2 , m 2 D o , m 2 + ) + A 3 / 2 {m 2 D .+ , m 2 D+ , m 2 )] /487rm^.+ 

(e(-3m; 1/2 + m-JJ 2 + 4m^)/(6 x 0.1361)) 2 (m^ + - m 2 D+ f 
A 3 / 2 (m 2 D *+ , m 2 D0 , m 2 + ) + A 3 / 2 (m 2 D *+ , m 2 D+ , m 2 
w (9.48 ± 0.43) x 10" 5 , 

where in particular we have used the supermultiplet prediction g VVP /g VPP = 2/m D *+ 
and we have confidence in our prediction to this order. This finding implies that the 
dominant decay modes in D* + decay are the PP channels we derived and as such 
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approximate well to the full width. Thus the radiative branching fraction for D* + — > 
Z} + 7 is relatively small. This is not inconsistent with recent measurements by CLEO 
which measured a fraction of (1.1 ± 1.4 ± 1.6)% [Q. However, our prediction does 



conflict with other theoretical models fL3| , |24| , pq , p6| . Branching fraction calculations 
for the two PP channels yield 

r D ^ D ^+/r D * + ^ n « 68.8% 

rD*+^D+7ro/r_D*+^all ~ 31.2% 

which compare well with other models and the experimentally determined results 
from CLEO (1992): 

AV+^ a ii = (68.0 ±1.4 ±2.4)% 

/r^+^aii = (31.0 ±0.4 ±1.6)%. 

We are able to employ similar methods in the decays of the D*° vector meson. 
In this instance, the possible PP decay channels are restricted by phase space to 
D*° — > D°7i°. We calculate this width to be 27.2±0.7 KeV, where we used go^D ^ = 
(0.1537 ± 0.002) (2007.1) 1 / 2 / v / 2. To determine the radiative width D*° D°j we 
apply relation ( ^9|) along with a small correction for the change in phase space to 
derive 

r D * ^ D o 7 « 3672 x r D *+^ D+7 = (20.1 ± 1.1) KeV. 

Thus the total D*° width is (47.3 ± 1.3) KeV, although the result is sensitive to the 
supermultiplet prediction g vvp (D*) / g VPP (D*) = 2/m D *. Consequently, we predict 
the following branching fractions 

r D *o^o 7r o/r D ,o_ all w (57.5 ±2.2)% 
r D *o_> D o 7 /T 0*0^x1 ps (42.5 ±2.6)% 

which are in fair agreement with the CLEO data: 

r D ,o_ D o w o/r D ,o^ all = (64 ±2.4 ±4.5)% 
r D .o_ D o 7 /r D »o_^ all « (36 ±2.4 ±4.5)%. 

6 Conclusions 

This study of two-body meson decays has shown that the supermultiplet method 
unifies meson decays quite well, even for the light quarks. The most significant finding 
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is that the coupling between mesons is susceptible to symmetry breaking mechanisms, 
but in a regular way, allowing us to successfully extrapolate to decay rates for other 
processes. In particular, the methods are readily applicable to heavy quark examples, 
as highlighted by our examination of D* processes. 
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Table captions. 

Table 1 Results for g VPP determination. 
Table 2 Results for g vvp determination. 



Figure captions. 

Figure 1 Quark line or duality diagram. 
Figure 2 Vector meson dominance. 

Figure 3 Support for the supermultiplet symmetry condition 2g vpp = mig vvp in 
the decays — > K°K° and — > pn. 
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